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1. Introduction 



The variety of accessible final states in B decays to two mesons provides an abundant source 
of information on CP violation and flavour-changing processes. When the final state consists 
of two vector mesons, an angular analysis of the vector mesons' decay products also provides 
insight into the spin structure of the flavour-changing interaction. For the V — A coupling of 
the Standard Model, a specific pattern of the three helicity amplitudes is expected [1], such 
that the longitudinal polarisation fraction fi should be close to 1. Since fi ~ 0.5 was first 
observed [2, 3] for penguin-dominated strangeness-changing decays, many theoretical papers 
addressed the question whether this result could be explained as a strong- interaction effect, or 
whether it could be reproduced within specific "New Physics" scenarios [4, 5, 6, 7, 8, 9, 10, 11, 
12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22]. 

In this paper we revisit this question using the QCD factorisation framework [23, 24] to 
deal with the strong interaction in the amplitude calculation. Our study goes beyond previous 
ones in several respects. On the theoretical side we provide the first complete results for the 
hard-scattering kernels relevant to vector- vector (VV) final states, correcting several errors in 
the literature. (In fact, the only correct calculation is [6].) We also provide a more detailed 
discussion of the factorisation structure and power counting for the various amplitudes. It 
seems to have escaped attention so far that, contrary to the longitudinal polarisation amplitude 
and those relevant to PP and PV final states, the transverse polarisation amplitudes do 
not factorise even at leading power in the heavy-quark expansion. This, together with the 
high sensitvity to penguin weak-annihilation [6], implies that the calculation of polarisation 
observables stands on a much less solid footing than the calculation of B — > PP, PV decays. 
On the phenomenological side, we provide estimates for all B — > VV decays (including B s 
decays) and for all parameters that enter the angular analysis. Previous studies concentrated 
on single or a few decay modes and considered the longitudinal polarisation fraction Jl only, 
making it often difficult to distinguish general patterns from the consequences of particular 
parameters choices. 

The organisation of the paper is as follows: In Section 2 we summarise the definitions for 
the helicity amplitudes, angular variables and polarisation observables. The calculation of 
the B — > VV decay amplitudes in the QCD factorisation framework is briefly reviewed in 
Section 3. We then discuss a few aspects of the transverse polarisation amplitudes that allow 
for an understanding of the main characteristics of B — > VV phenomenology. One important 
conclusion from this discussion is that the analysis of B — > VV decays will be much less 
rigorous and much more uncertain than the corresponding analysis of B — > PP and B — > VV 
modes [25]. The technical results of the calculation are summarised in an Appendix. Section 4 
provides the list of input parameters, an overview of the flavour amplitude parameters with 
theoretical uncertainties, and a classification of the 34 B — > VV decay channels, which guides 
the subsequent numerical analysis. We begin the analysis in Section 5 with a discussion of 
branching fractions, CP asymmetries and polarisation observables of the nine tree-dominated 
decays. Among these the four colour-allowed modes can be well predicted. In particular, we 
show that the time-dependent CP asymmetry measurement in B° — > p + p~ leads to one of the 
most accurate determinations of the CKM angle 7. In Section 6 we turn to the 14 colour- 
allowed penguin-dominated decay modes. It will be seen that theoretical calculations allow 
for large transverse polarisation within large uncertainties. This suggests to determine the 
transverse penguin amplitude from data using the well-measured <\>K* modes. This approach 
is used to sharpen the predictions for the remaining decay modes in this class. The analysis 
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concludes in Section 7 with a brief discussion of the remaining penguin-dominated modes, and 
decays that occur only through weak annihilation. Section 8 summarises our main results and 
conclusions. 



2. Helicity amplitudes and polarisation observables 

We consider a B meson with four- momentum ps and mass tub decaying into two light vector 
mesons Vi(pi,r/*), V2(p2>0 with masses mi^ of order Aqcd- The decay amplitude can be 
decomposed into three scalar amplitudes Si,2,3 according to 

Ab^vm = ir)*^*» f Si 9lw - S 2 ™^ + S 3 ie^V^V) . (1) 
V rn B VrV2j 

with convention £0123 = 1- Alternatively, one can choose a basis of amplitudes describing decays 
to final state particles with definite helicity 

im 2 B f S 2 



A± = A(B - Vi(pi,i4)V 2 (p2,4)) — i (Si T S 3 ). 

or use the transversity amplitudes, where .A± are replaced by A\\ = {A+ + A-)/y/2 and 
-4_i_ = (.4+ — A-)/V2, corresponding to linearly polarised final states. We choose p 2 to be 
directed in positive z-direction in the B meson rest frame, and the polarisation four-vectors of 
the light vector mesons such that in a frame where both light mesons have large momentum 
along the z-axis, they are given by e± = i]^ = (0, ±1, i, 0)/\/2, and 6q = p 2 l /m2, r/^ = p^/mi. 
Here and thoughout the paper we neglect corrections mj 2 /wi| quadratic in the light meson 
masses. Thus p^ 2 = mBn^/2 with n± = (1, 0, 0, ±1). Our conventions imply the identity 

i£ tJl vpae±n+ = (Tl) {n +fi e* ±u - n +l/ e* ±ll ) . (3) 

Experimentally, the magnitudes and relative phases of the various amplitudes are extracted 
from the angular distributions of the vector resonance decay products. The full angular de- 
pendence of the cascade where both vector mesons decay into pseudoscalar particles is given 

by [i] 



dT 



OC I Aq 1 2 COS 2 0i COS 2 #2 + J Sin2 01 Sin2 $2 (\A+\ 2 + \A-\ 2 ^j 



cfcos Videos "&2dip 

-cosi?isini?icosi? 2 sini? 2 [Re (e~ ilp A A* + ) + Re (e +ilp A A*_)] (4) 
+ i sin 2 i?i sin 2 2 Re {e 2i ^A+A*_) , 

where ip measures the angle between the decay planes of the two vector mesons in the B meson 
rest frame, and 0^2 are the angles between the direction of motion of one of the Vi 2 - ► PP 
pseudoscalar final states and the inverse direction of motion of the B meson as measured in 
the Vi t 2 rest frame, see Figure 1. The omitted proportionality factor is such that one obtains 
the decay rate T(B — ► V1V2) when cos$i and cos $2 are integrated from —1 to 1, and ip from 
to 2tt. 
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Figure 1: Decay kinematics in the rest frame of V\. 



Thus, any given B — > VV decay allows us to define five observables corresponding to the 
three magnitudes and two relative phases of the helicity amplitudes, or the five angular co- 
efficients in (4). In experimental analyses, observables are preferably defined in terms of the 
transversity amplitudes as they have definite CP transformation properties. A typical set 
of observables consists of the branching fraction, two out of the three polarisation fractions 
fL,f\\,f±, and two phases cj>±, where 

^ = ■„ .2 'mi',, ,2 ' C = ar ^- ( 5 ) 

A n r A- L4.i -4- \ A , r 11 An 





A o,\u 


2 


\A \ 2 + 


A \\ 


2 


+ l-4±| 2 



It is conventional to combine the five observables of some B — > VV decay with those of its 
CP-conjugate B decay, and to quote the ten resulting observables as CP-averages and CP- 
asymmetries. We denote B decay helicity amplitudes as Ah and define the corresponding 
transversity amplitudes as Ann = (A- ± A+) / V2, so that A\\/± = A\\/± in the absence of CP 
violation. Observables f B ,4>h are then defined as in (5), and CP averages and asymmetries 
are calculated by 

fh = Uti + tf), ^P = f^k (6) 



(7) 



z 4+4 

(h = L, [|, _L) for the polarisation fractions and 



— J.B 



B - A(j) h (mod 2vr) 

h + A<p h (mod 2tt), -| < A(f> h < | 



(/i = ||, -L) for the phase observables (fit and A^. The implicit definition (7) ensures that the 
CP-averaged phase is the geometrical bisection of the acute angle enclosed by <p B and (fi B ; the 
magnitude of this angle is 2 \A(f>\. More explicitly, the averaged quantities can be obtained as 

1 

2 
1 



TT 



A^ h = -(<P B -cp B )+7r-e( <ft-ffi -tt) (mod 2vr). 



(8) 



Our phase convention for the amplitudes and definition of observables is compatible with 
that used in the relevant publications of the BaBar and Belle collaborations (for example 
in [26, 28, 29, 30]), except for the sign of .Ao relative to the transverse amplitudes, which leads 
to an offset of ir for <f>\\/±- We favour the above convention, because it implies (f>\\ = 4>± = 
and A(p\\ = A<p± = at leading order, where all strong phases are zero. 
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3. B — ► VV amplitudes 

The decay amplitudes follow from the matrix elements {ViV2\H c ff\B) of the effective Hamilto- 
nian (conventions as in [31]) 



* p=u,c 



with D G {d, s} and = V p bV* D . A quark model [1] or naive factorisation analysis indicates 
a hierarchy of helicity amplitudes 



>+h.c. (9) 



i=3,...10,7 7 ,8g 



m 6 



. MqcdV 
\ rn b J 



for S meson decays. (For .B decays exchange — <->+.) This is a consequence of the left- 
handedness of the weak interaction and the fact that high-energy QCD interactions conserve 
helicity 

In naive factorisation one considers only the four-quark operators in 7i c R and approximates 
their matrix elements by the matrix elements of two currents [32]. The helicity amplitudes 
y are proportional to 

A h VlV2 = ^§ (V^\(q s b) v ^ A \B qs )(V 2 h \(q q ') v \0) (11) 

in this approximation. Evaluating this expression (conventions for the form factors as in [33]) 
we obtain 



A° VlV2 = l -^ m 2 B f V2 A^ Vl (0), A± iV2 = l -^ mB m 2 f V2 Fi^(0) (12) 



with the definitions 

*£- v V) - (1 + Z±)A*-**tf) T (1 - ?±)V™tf). (13) 
me me 

The transverse amplitudes A Vi y 2 are suppressed by a factor m^jvnB relative to A Vi y 2 . In 
addition, the axial-vector and vector contributions to F^^ Vl (0) cancel in the heavy-quark 
limit, due to an exact form factor relation [33, 34]. Thus F_/Aq ~ 1, F + /Ao ~ C(Aqcd/^_b), 
and (10) follows. 

The dominance of the longitudinal amplitude indicated by (10) leads to the well-known 
expectation that fi should be close to unity. Experimental data for penguin-dominated B 
decays is in conflict with this expectation thus motivating theoretical studies beyond the naive- 
factorisation approximation. 

3.1. The QCD factorisation approach for B — > VV 

We use the QCD factorisation approach [23, 24] to compute the matrix elements (V\V2\Qi\B) 
of the effective Hamiltonian. In this framework they can be expressed (at leading power 
in an expansion of the amplitude in Aqcd/^b) in terms of form factors, meson light-cone 
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distribution amplitudes and perturbatively calculable hard scattering kernels. In condensed 
notation, the factorisation formula reads 



{ViV 2 \Qi\B) = (F B ^ Tl * f V2 $ V2 + [Vi <- V 2 \) + T! 1 * f B >5> B * f Vl $ Vl * Jy 2 $v 2 , 



(14) 



where the star products imply an integration over light-cone momentum fractions. In addition 
the framework contains estimates of some power corrections, which usually cannot be computed 
rigorously. 

We follow closely the scheme developed in [25] for B — ► PP, PV decays to match contri- 
butions to the hard-scattering kernels T?' 11 on terms involving products of flavour coefficients 
o%' (V1V2) and factorised matrix elements Ay iV2 . The longitudinal amplitude h = can be 
deduced from the results given in [25]. For the analysis of the present paper we calculated 
the transverse helicity amplitudes. In the following we describe the basic results and main 
differences with respect to the longitudinal amplitude; the expressions for the hard-scattering 
functions are given in the Appendix. 

Non-leptonic decay amplitudes are sums of products of CKM factors, Wilson coefficients 
from (9) and matrix elements (14) of operators with different flavours. It is convenient to 
organise the amplitudes according to flavour. Thus, one writes for example, 



V2A h B -^ p0K *- = £ A« U h pR , [5 pu (of + + of* + <£w + Pi* + 01 



,p,h 
3,EW 



p=u,c 



-J- A h - 



V«2 + ^;ew }■ (is) 



In naive factorisation, the flavour coefficients a?' h (V\V2) are linear combinations of Wilson 
coefficients Cj. In QCD factorisation, they include non-factorisable loop effects and spectator- 
scattering. The (3f ,h coefficients parameterise weak annihilation amplitudes. The decompo- 
sition of the amplitudes for the 34 VV final states in terms of these quantities follows from 
the PV expressions given in [25] with obvious replacements of pseudoscalar by vector mesons. 
The ai relate to the coefficients used in the older factorisation literature as follows (helicity 
indices and V1V2 arguments suppressed): 



\ = a\- r£a\ 



a 2 =a 2 «3,EW = a 9 + a 7 «4,EW = a 10 ~ r x °8> 



where we have used the notation 

rl = . (17) 

m b f v 



The explicit expressions for the negative-helicity coefficients a?~ and the transverse weak 
annihilation amplitudes are collected in the Appendix. Beyond leading order the a?~ are sums 
of vertex corrections, penguin contractions, and spectator-scattering contributions, see (53). 
Most of the relevant hard-scattering functions have been calculated before. In [4, 5, 12, 13, 16] 
results for all three of these contributions have been given. We do not find agreement with 
these results, however. As far as we understand, the origin of the discrepancy is that the 
authors of these papers use an incorrect projection on the light-cone distribution amplitudes 
of transversely polarised vector mesons, which neglects the transverse-momentum derivative 
terms in (52). An exception is [13], which does state the correct projector, but the results 
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still differ from ours, particularly for the spectator-scattering contributions. Kagan [6] has 
calculated the QCD penguin contractions, spectator-scattering terms as well as the weak an- 
nihilation amplitudes, but did not consider the vertex contractions. We confirm his results on 
the penguin contractions and QCD penguin annihilation, for which explicit expressions were 
given in the paper. 

3.2. Anatomy of transverse amplitudes 

The NLO calculation of the negative-helicity amplitude is quite similar to the calculation of 
the longitudinal amplitude. The result exhibits, however, some qualitative differences which 
have important consequences for the phenomenology of B — ► VV decays. In this section 
we explain the non-factorisation of the negative-helicity amplitude; that the positive-helicity 
amplitude cannot be calculated in an analogous manner; that the amplitude hierarchy (10) is 
violated by electromagnetic effects; that penguin annihilation is comparatively more significant 
for transverse polarisation than longitudinal polarisation penguin amplitudes. 

3.2.1. Non-factorisation of spectator scattering 

The factorisation formula (14) contains two structurally different terms, the first of which is 
dominated by soft interactions within the B — > M\ transitions. These are absorbed into the 
QCD form factor. The second term stands for interactions where a hard (more precisely, hard- 
collinear) interaction with the spectator quark in the B meson takes place. Both terms are of 
the same order in the heavy-quark expansion. This remains true for the transverse polarisation 
amplitudes, but now one finds that the convolution integrals over the light-cone distribution 
amplitudes are logarithmically divergent due to the occurence of the integral 



see (61), (62). This endpoint divergence at x = 1 signals that the presumed factorisation 
of spectator-scattering does not hold even at leading power. A similar effect occurs in the 
B — > PP,PV and longitudinal B — > VV amplitude only when one attempts to calculate 
power corrections by applying the light-cone projection including twist-3 terms. It is perhaps 
not surprising that this divergence is obtained at leading-power for the transverse amplitude, 
since the entire amplitude is formally a twist-3 term. (This is the origin of the power suppression 
of the transverse amplitudes relative to the longitudinal amplitude.) Factorisation-violation 
at the leading power implies that the calculation of transverse polarisation amplitudes is on 
a much less solid footing than of the other amplitudes, and often should be considered more 
as an estimate. In practice, we find that the non-factorisation of spectator-scattering is only 
significant for the colour-suppressed tree and the flavour-singlet QCD penguin amplitudes, 
where the (regulated) divergent integral is multiplied by a large Wilson coefficient. The fact 
that there are endpoint divergences in the spectator-scattering contribution to the transverse 
helicity amplitudes has been observed in previous calculations, but its significance for the 
theoretical status of the factorisation approach and its phenomenological implications have 
not been sufficiently emphasized. 




(18) 
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3.2.2. The positive-helicity amplitude 

The calculation of the kernels T/ in (14) can be interpreted as matching the operators Q% to 
four-quark operators with field content [xx][£^u] i n soft-collinear effective theory [35]. The field 
X describes collinear quarks moving in the direction of V2, the meson that does not pick up the 
spectator quark, and satisfies fi+x = 0. The leading quark bilinears that have non-vanishing 
overlap with (V2I are 

X^-(lT75)x, xj4-7±(l±T5)x- ( 19 ) 

The subscript _L denotes projection of a Lorentz vector on the plane transverse to the two 
light-cone vectors n=p. The first operator overlaps only with the longitudinal polarisation state 
of V2, the second only with a transverse vector meson. However, the second operator is not 
generated by the V — A interactions of the Standard Model, at least up to the one-loop level. 
Hence the transverse amplitudes arise from power-suppressed operators 0± = [X-D±x][£^]- 
(In terms of the light-cone projector (52) this statement implies that the leading term in the 
first line does not contribute for V — A interactions, leaving the twist-3 terms in the second and 
third line. Since 0± contains transverse momentum derivatives, one must keep the transverse 
momenta of partons collinear to V2; this explains why the transverse- momentum derivative 
terms in the projector are required.) The left-handedness of the weak interactions implies that 
operators of this form contribute only to the negative-helicity amplitude. The positive-helicity 
amplitude appears first in yet higher-dimensional operators such as [xD±x][£D±h v ]. To match 
to such operators one must keep the transverse momentum of the quark lines collinear to both 
mesons non-zero. Such a calculation has not yet been done, and therefore all calculations of 
the positive-helicity amplitude in the literature must be regarded as incomplete. It is even 
possible that for the positive-helicity amplitude no useful factorisation formula holds even for 
the non-spectator-scattering terms in (14). 

It follows that the positive-helicity amplitude is power-suppressed relative to the negative 
one, and should be set to zero in the absence of any consistent calculation of this power 
correction. Within this approximation, A+ = A- = 0, there are only two rather than four 
independent polarisation observables, since 

f\\ = f±, <f>\\ = <f>±- (20) 

Similarly identities for the corresponding CP asymmetries hold. It should be noted that these 
identities are non-trivial consequences of the V — A nature of the weak interactions and of 
factorisation, and it is therefore worthwhile to test them experimentally. 

In our analysis we proceed as follows: we assume the naive-factorisation expression for the 
positive-helicity amplitude A+ and allow the form factor to vary within the range F^^ Vl = 
± 0.06. Thus, we allow a small variation of A+ around zero to estimate the error from 
neglecting this power correction. We note that QCD sum rule results for the form factors give 
F + values consistent with zero [36, 37]. 

3.2.3. Violation of the amplitude hierarchy 

In the previous paragraphs we explained the origin of the amplitude hierarchy (10). How- 
ever, when electromagnetic effects are included, a transverse polarisation amplitude can be 
generated by a short-distance transition to a vector meson and a photon with small virtuality 
which subsequently converts to a vector meson [20]. This transition is enhanced by a factor 
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(m^/Agco) 2 due to the large photon propagator, resulting in the parametric relation 

Ao :A- : A+ = 1 : : a em . (21) 

A QCD 

Thus, formally, the negative-helicity amplitude is leading in the heavy-quark limit. Techni- 
cally, this is related to the existence of the operator [yl 7± ][^/i„], which contributes only to the 
transverse amplitude, and which is enhanced relative to [xx][£^-u]- In the Standard Model the 
leading effect of this type involves the electromagnetic-dipole operator in the effective Hamilto- 
nian. As a consequence the colour-allowed electroweak penguin amplitude a-j^w * s completely 
different from its naive-factorisation value. Our calculations include this contribution, which 
has already been discussed specifically in [20]. 



3.2.4. Penguin weak annihilation 

Weak annihilation is a power correction not included in (14), since it does not factorise due to 
endpoint divergences in the convolution integrals. The effect is often estimated by a parame- 
terisation suggested in [31], where the endpoint-divergences are regulated by a cut-off. In this 
model one finds that the most important annihilation effect is a penguin annihilation ampli- 
tude that is phenomenologically indistinguishable from the QCD penguin amplitude. These 
general observations also hold for the transverse polarisation amplitudes. In particular, the 
weak annihilation contribution to the negative-helicity amplitude is a power correction relative 
to the leading, factorisable contributions to this amplitude. Yet, as found in [6], the effect is 
numerically much larger than in B — > PP, PV decays, and perhaps so large that a theoretical 
calculation of the negative-helicity QCD penguin amplitude is no longer possible. 
To explain this point, we consider the QCD penguin amplitude 



r>h _ Ah 



<4 + /#l, (22) 



and compare the h = and h = — amplitudes. Here a\ are the QCD penguin contributions, 
and /?3 the penguin annihilation contributions. For the longitudinal amplitude /3§ is Aqcd/^-6 
suppressed relative to aq, but it turns out that numerically the largest effect arises from a 
(Aqcd/^6) 2 term, which has a large colour factor and Wilson coefficient. This particular con- 
tribution is not suppressed by the extra factor of Aqcd/w-;) in the negative-helicity amplitude. 
One finds A Vi y 2 (3^ « Ay V /3%, while Ay V a^ <C Ay y a° due to the power suppression of 
A Vi y 2 with respect to A Vi y 2 . Thus, relative to a\ the numerical effect of (3% is a factor of 
mb/AQCD larger for the negative-helicity amplitude than for the longitudinal amplitude (but 
still power-suppressed since the suppression was (Aqcd/^6) 2 for h = 0). 

The following numerical estimates illustrate this point. We consider the p = c, h = 0, — 
helicity amplitudes for pK*, and also the ttK amplitude for comparison. The imaginary parts 
of the amplitudes are neglected, since they are not important for this discussion. We then find 

al^K) + fc{*K) = -0.09- {0.02 [-0.01, 0.05]}, 

af(pK*) + l3f( P K*) = -0.03- {0.00 [-0.00, 0.00]}, (23) 

a c 4 -(pK*)+ !3 c f{pK*) = -0.05- {0.03 [-0.04,0.10]}. 

The numbers in curly brackets refer to (3^ ,h . The first number in brackets is the default 
value, while the interval provides the range allowed by the parameterisation adopted in [31]. 
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We observe a (presumably accidental) cancellation in the longitudinal annihilation amplitude. 
What is significant is the difference in the range of the interval relative to a c £ h for the negative- 
helicty VV amplitude vs. the PP amplitude. In particular, the annihilation contribution fit 
may be significantly larger than the QCD penguin amplitude at for the negative-helicity 
amplitude. We can also compare the h = — and h = amplitudes, 

P~ A pK* at + fit 0.05 +[-0.04, 0.10] 

P* * 3£ af * 0^2 ' (24) 

where we used A~ Km /A^ Kt « 1/4. This shows that P~ could be as large as P°, if annihilation is 
maximal. Thus, for penguin-dominated decays, a longitudinal polarisation fraction fi around 
0.5 is not ruled out. 

Let us summarise these and a few further observations on the role of weak annihilation in 
B — > decays: 

1) The annihilation contribution to the longitudinal penguin amplitude is small, perhaps 
due to an accidental cancellation. 

2) The annihilation contribution to the negative-helicity penguin amplitude can (but need 
not) be very large, possibly leading to significant transverse polarisation in penguin- 
dominated decays. 

3) No such enhancement is observed for the annihilation contribution to the tree amplitudes, 
hence tree-dominated decays should be predominantly longitudinally polarised. 

4) We also calculated the weak annihilation contribution to the positive-helicity amplitude, 
and find that the large contribution to fi% is absent. Hence there is no evidence for large 
corrections to (20) even for penguin-dominated decays. 

It should be clear that these statements assume that the parameterisation adopted in [31] 
reproduces correctly the qualitative features of the weak annihilation amplitudes. 

4. Input and overview 
4.1. Input parameters 

The values of the Standard Model and hadronic input parameters are listed in Table 1. When 
we compare B — ► VV modes to decays with pions in the final state, the additional pion 
parameters are / w = 131 MeV, a 2 (2GeV) = 0.2 ± 0.15, ff n {0) = 0.25 ± 0.05. (The light 
quark mass values reported in the Table are only needed for the computation of the pion 
decay amplitudes.) Relative to the analysis of PV final states [25] we have implemented 
several minor parameter modifications (Wolfenstein parameter A, \V c b\, B meson lifetimes 
and decay constants, Gegenbauer moments of light-meson light-cone distribution amplitudes), 
which reflect new measurements or improved calculations, but individually have little impact 
on the calculation of non-leptonic decay amplitudes. A more important change concerns the 
treatment of \V u b\ and the B meson parameter A#, where we (roughly) stick to the same 
ranges as before, but choose smaller default values. These values lead to a good agreement 
of theoretical calculations with the observed B — > tttt transitions as already noted in [25]. 
Regarding the value of \ V u t>\ we note that our default value correspnds to the one that is favoured 
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QCD scale and running quark 
Aqcd = Ai± m b (m h ) 


masses [GeV] 
m c 


m s (2GeV) 


m q /m s 


0.225 4.2 


1.3 ±0.2 


0.09 ±0.02 


0.0413 


CKM parameters 

A \V cb \ 


1 v ub /v cb | 


7 




0.225 0.0415 ± 0.0010 0.085±g;gf| 


(70 ± 20)° 




B meson parameters 
lifetime r [ps] 
decay constant fs [MeV] 
X B [MeV] 


B- B° 
1.64 1.53 

210 ±20 

200^0 


B s 
1.46 
240 ± 20 


Light meson decay constants and Gegenbauer moments 

p K* LO (f> 


f [MeV] 209 ± 1 
f ± [MeV] 150 ± 25 
ai, ai t ± 
«2, a 2 ,± 0.1 ±0.2 


218 ±4 
175 ± 25 
0.06 ± 0.06 
0.1 ±0.2 


187 ±3 
150 ± 25 


0.0 ±0.3 


221 ±3 
175 ± 25 


0.0 ±0.3 


Form factors for vector mesons at q 2 = 
B -> p B -> K* B ^lo 


B S ^K* 


B s ^4> 


Aq 0.30+^;^ 0.39 ± 0.06 
F_ 0.55 ± 0.06 0.68 ± 0.07 
F + 0.00 ± 0.06 0.00 ± 0.06 


95 +0 - 07 
0.50 ± 0.05 
0.00 ± 0.06 


0.33 ± 0.05 0.38lg;J§ 
0.53 ± 0.05 0.65ig;^ 
0.00 ± 0.06 0.00 ± 0.06 



Table 1: Summary of theoretical input parameters. All scale-dependent quantities refer to 
ji = 2 GeV unless indicated otherwise. m q = (m u + m&)j1. 

by exclusive semi-leptonic b — > u transitions, which is smaller than the one from the inclusive 
decays (see [38] for the most recent discussion). We shall see below that the tree-dominated 
B — > pp modes also support this small value unless the B — ► p form factors are unacceptably 
small, so all exclusive decays (semi-leptonic, non-leptonic PP, PV and VV modes) seem to 
consistently favour small | V ub \ . The default value we adopt for A# is significantly smaller than 
the value obtained from QCD sum rule calculations [39, 40], which we use to define the upper 
limit of this parameter's range. Some of the longitudinal and negative helicity B — > V form 
factors have changed considerably since the publication of [25] due to the update of the QCD 
sum rule calculation [41]. Our new values follow [41], but in some cases a smaller form factor is 
adopted to improve the description of data. The smaller values are compatible with [41] within 
theoretical errors; in these cases, however, the theoretically allowed parameter range becomes 
asymmetric around the default value. The positive-helicity form factors are set to 0.00 ± 0.06 
(see Sect. 3. 2. 2). The renormalization scales are treated as in [25] and p is varied from m b /2 to 
2m b . The Wilson coefficients Cj are tabulated in [31]. 

In addition to the well-defined hadronic parameters, the predictions of QCD factorisation 
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Parameter 



h = 



ai(pp) 

«2(PP) 

Pi(pp) 

Upp) 


Q4+ - 10 

u - Ji -0.25 

n 05 +0 - 03 

u - uo -0.06 

-n n2 +0 - 02 

u - uz -o.oi 


+ 
+ 
+ 
+ 


( o.o2±°:8|)i 
(-costal')* 

(-0.02±g : oi)i 




1 14+0.30 , 

-o.i9i8:? 6 7 + 
o.oii8:81 + 
o.ooi8:88 + 


( o.o4!8:D; 
;-o.i7l8:|)i 
: o.ooi8:81), 
: o.ooi8:88), 




a%(K*cf>) 
c%(K*<t>) 

0l{pp) 


u.uuo_ 004 
-0 024+ - 004 

-n rm +a007 

u.uoo_ 007 

-o ono +0 - 002 

U.UUU_q QQ2 

-o.oo5±g : g82 


+ 
+ 
+ 
+ 
+ 


(-0.00118:88^ 
(-o.oi5i8:88 5 

on+0-005 

I u - uli -0.004 

f 001+ 0001 

v u - uui -0.002 

f nnn2 +0 - 003 

l u - uuz -0.006 


H 
H 
)i 
H 
)i 


-o.oo5i8:811 + 
-o.o49i8:81^ + 

-o 047+ 0017 + 

U - U4 ' -0.017 ^ 
u.uoi_ 065 -f 

u.uui_ 001 -+■ 


;-o.ooii8:81^ 
(-o.oi7i8:81^ 
(-o.oo2i8:81^ 
( o.o3ii8:^ 
( o.oooi8:881 


)i 
)i 
)i 
)i 
)i 


<4,EW(K*P) 
<Ew(*"» 


-o.oo8i8;88i 
-0.00218:88^ 


+ 
+ 


(-0.00018:888 
( o.ooii8:881 


)i 
H 


o.oi5i8:88^ + 
o.oo2i8:88? + 


(-o.oo3i8:88^ 
( o.ooii8:88? 


)i 
)i 



Table 2: Overview of longitudinal and negative-helicity amplitude parameters. 



depend on the model parameters Xh and Xa, Xl (see [31] and the Appendix for their defi- 
nition). In contrast to all previous QCD factorisation calculations there is model-dependence 
even at leading power in the heavy quark expansion for the transverse amplitudes due to an 
endpoint divergence in spectator scattering (see Sect. 3.2.1). This is parameterised by 

X H = (l + Q H e i ^)\n r ^-- A^ = 0.5GeV, (25) 

with qh = by default, a range defined by qh < 1, and an arbitrary phase ipu- Xa related 
to weak annihilation is defined in the same way, but here we use qa = 0.6 e~* 40 by default. 
The motivation for this choice will be explained in the context of penguin-dominated decays. 
In B — > VV decays the parameter Xa is only relevant for the negative-helicity penguin am- 
plitudes due to the near-cancellation of longitudinal weak annihilation (see Sect. 3.2.4). For 
completeness we note that Xl is evaluated by an equation similar to (25) with ln(m,B / A-h) 
replaced by ms/A^, but Xl is never numerically relevant in the amplitude calculation. 

In the end all parameter (Standard Model, hadronic, model) uncertainties are added in 
quadrature, except for the CKM parameters, which are separated, because the dependence 
on the CKM parameters \V u b\ and 7 is interesting: if larger than the hadronic error for some 
observable, this observable may be useful to determine \ V u b\ or 7. In general, the first "error" on 
a quantity will provide the dependence on CKM parameters; the second gives the "theoretical 
uncertainty" . 
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4.2. Flavour amplitudes 

The helicity decay amplitudes such as the example (15) are composed of CKM factors, the 
factorisable coefficients (12) and the "flavour amplitudes" af, (5^. The flavour amplitudes 
correspond to the colour-allowed (colour-suppressed) tree amplitude (a^), the QCD penguin 
amplitudes (a^' ), the QCD singlet-penguin amplitudes (a^ ,h ), and the colour-allowed (colour- 
suppressed) electroweak penguin amplitudes a^ w («4E W )- There also exist tree annihilation 

Pi 2 , penguin annihilation (5^\ and electroweak penguin-annihilation amplitudes /?f 4ew The 
result of our calculation of these amplitudes (except for the irrelevant electroweak annihilation 
amplitudes) is summarised in Table 2. Most of the later analysis of branching fractions, 
CP asymmetries and polarisation observables can be reproduced by inserting these numerical 
estimates into the expressions for the decay amplitudes in terms of flavour parameters in 
the appendix of [25]. The flavour parameters depend on the final and initial state, but this 
dependence is rather small and may be ignored for rough estimates. An exception is the 
negative-helicity electroweak-penguin amplitude, since the power-enhanced electromagnetic 
contribution depends quadratically on the light-meson mass. The Table gives the numbers 
for the pK* final states, where the electroweak penguin amplitudes have the most significant 
effects [20]. 

Let us point out the most important features of the B — > VV amplitudes related to the 
general discussion in the previous section. The longitudinal QCD penguin amplitude a^'° is 
rather small, similar to the VP or PV penguin amplitudes. However, for VV the QCD pen- 
guin annihilation amplitude f3^'° is strongly suppressed, and irrelevant, in marked difference to 
the case of B — > PV decays. A striking result for the negative-helicity amplitudes is the value 
and large uncertainty of the colour-suppressed tree amplitude, and to some extent even of the 
colour-allowed tree amplitude. This reflects the non-factorisation of spectator-scattering. The 
same effect is also responsible for a larger uncertainty and, possibly, large value of the QCD 
singlet-penguin amplitude, which may therefore be relevant to the (f>K* modes. As has already 
been discussed in some detail, the QCD penguin annihilation amplitude is large, perhaps 
larger than af^~ , which is evident from the Table. Finally we note that the negative-helicity 
electroweak penguin amplitude o^'ew has a different sign from the corresponding longitu- 
dinal amplitudes, which is a consequence of the additional, power-enhanced electromagnetic 
contribution. 

4.3. Classification of decay modes 

We conclude this overview section with a classification of the total of 34 B~, B° and B s decay 
channels into two light vector mesons according to the reliability of the calculation of various 
observables. This classification is motivated by the observation that the peculiarities of the 
transverse-helicity amplitude calculation analysed in Section 3.2 severely limit the reliability 
of QCD factorisation for many observables. In effect, of all the transverse amplitudes, only the 
colour-allowed tree and electroweak penguin amplitudes can be calculated with some accuracy. 

• Colour-allowed tree- dominated AD = 1 decays. Most observables are amenable to cal- 
culation, the exception being CP asymmetries of polarisation observables, which always 
involve transverse QCD penguin amplitudes. Only four decays belong to this class, 
namely B° — ► p + p~ , B~ — > p~p°, B~ — > p~u>, and B s — > p~ K* + . 

• Colour- suppressed tree- dominated AD = 1 decays. The non-factorisation of spectator 
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scattering in the transverse amplitude precludes a reliable calculation of polarisation 
observables for these decays. However, if the longitudinal amplitude is still dominant, 
predictions for the CP-averaged branching fractions and Acp can be obtained, and the 
longitudinal polarisation fraction is expected to be close to 1. The five modes B° — > p°p°, 
B° -► pV B° -> coco, B s -> p°K*°, and B s -► coK*° fall into this category. 

• Penguin- dominated decays. In these modes, no polarisation observables can be calculated 
reliably from theory alone because of penguin weak-annihilation effects. As transverse 
and longitudinal contributions cannot be excluded to be of similar magnitude, even 
branching fraction and CP asymmetry predictions will suffer large uncertainties. The 
eleven AS* = 1 modes B — ► pK* ,ooK* ,<pK* , B s — ► K*K*,<p<p with branching fractions 
in the upper 10 -6 range, and the three AD = 1 modes B° -► K*°K *, B -► K*°K*-, 
B s — > </>if *° with small branching fractions belong to this class. 

• Electroweak or QCD flavour- singlet penguin- dominated decays. These decays are ex- 
pected to have very small branching fractions. They are difficult to predict, if the QCD 
flavour-singlet penguin amplitude plays a role. The five decays B~ — ► p~4>, B° — > p ^, 
5° — > and i? s — > l? s — ► w0 belong to this class. The B s decays in this class may 
exhibit a significant, perhaps dominant, contribution from the doubly CKM-suppressed, 
colour-suppressed tree amplitude. 

• Pure weak annihilation decays. The six decays falling into this category, namely B° — > 
K*~K* + , B° — > f? s — > -B s — ► S s — > p°oo, and S s — > cjoj, are completely 
annihilation model-dependent, and only rough estimates of their branching fractions can 
be given. 

As we proceed with our analysis, we will discuss these categories in order for the tree-dominated 
and penguin-dominated decays. 

5. Tree- dominated decays 

The four tree-dominated colour-allowed modes are among the few B — ► VV decays that can 
be reliably calculated in QCD factorisation. They fully respect the helicity amplitude hierar- 
chy (10) and should exhibit predominantly longitudinal polarisation. Power corrections have 
limited impact, and the main sources of theoretical uncertainties are V u b and form factors. 
Penguin amplitudes are small, implying small direct CP asymmetries and the prospect of a 
precise determination of sin 2a from time-dependent B° — > p + p~ studies. On the contrary, 
the colour-suppressed tree-dominated decays have much smaller branching fractions, and the 
theoretical calculations are limited by large uncertainties in the colour-suppressed amplitude 
ct2, in particular in the transverse amplitude, where non-factorisation of transverse spectator 
scattering can spoil the predictions. In this section we quantify these expectations. 

5.1. Branching fractions and direct CP asymmetries 

We present the CP-averaged branching ratios and direct CP asymmetries in Table 3. 

It is interesting to note that the experimental data on the pp branching fractions exhibit a 
pattern similar to the corresponding -kit modes. The p + p~ and p~ p° modes have nearly equal 
branching fractions, and the p°p° has a larger branching fraction than naively expected. For 



13 



BrAv / 10~ 6 



Aqp I percent 



Theory Experiment Theory Experiment 

18.8±g-^±|-| (*) 18.2 ±3.0 0±g±g -8 ±13 

(*) 23.1+1;! -1±8±| +11 ±13 

O-S-as^ol 1.07 ±0.38 +28±|±|| n/a 

12.811;^ (*) lo.eit! -sliU +4±is 

0.2±^ ;[±g f < 1.5 no prediction n/a 
0-9±8:I±S:B <4-0 -29±g±» n/a 
25.2+|fif]' (*) n/a -3+^+j n/a 
1.5±S : g±?^ n/a ±27 + 5+3 4 7 n/a 
lf±8jgj Va -34+j°+g n/a 

Table 3: CP-averaged branching fractions and direct CP asymmetries of tree-dominated 5 — > 
decays. Experimental values are taken from [42, 43, 44, 45, 46, 47]. For numbers 
marked with an asterisk, the dependence on |V^,| and the form factor is obtained as 
described in the text. Errors are calculated as described in Section 4.1. 
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both, pions and p mesons, this is attributed to a larger colour-suppressed tree amplitude. In 
the factorisation framework this is realized if spectator scattering is the dominant dynamical 
mechanism behind the colour-suppressed tree amplitude. This favours the parameter choice 
adopted in Section 4 with small A# and small V u b and/or form factors. It is seen from the 
Table that the existing data is consistent with the theoretical calculation. 

The colour-suppressed decays are easily distinguished in the Table by their small branching 
fractions and large relative uncertainties. These uncertainties are dominated by the parameters 
for power corrections (mainly Xh relevant to spectator scattering) and there is little room for 
improvement from theory alone. The CKM error on the branching fractions is dominated 
by \Vvb\- Except for B° — > cop , the dependence on \ V u b\ can be extracted by assuming that 
the branching fraction is proportional to |V^| 2 . It is worth analysing the uncertainties of the 
colour- allowed decays in more detail, since they are dominated by V u b and form factors. Both 
can be expected to be known more accurately in the nearer future. Only the longitudinal form 
factors A^ Ml (0) are relevant here, since the transverse amplitudes contribute only a small 
amount to the branching fraction. To make this dependence explicit, we write 

2 

x (I8.8±8:*±i:|) • ID" 6 (26) 

for the p~p {) , and similarly for the p + p~ final state. For B~ — > uop~ and B s — ► K* + p~ 
we extract the default values of ^^"^(0) and Aq s ~* k (0), respectively. We then quote the 
number 18.8±g;|±|;| in Table 3, where the error includes all parameters but only the residual 
dependence on \V u b\ and the form factor. The bulk dependence can then be obtained by 
inserting into (26) whatever values of \ V u b\ and the form factor one prefers. We may conclude 
from the branching fractions of the colour-allowed decays that a significantly larger value of 
\V u b\ (than our default value \V u b/V c b\ = 0.085) is only compatible with data, if all form factors 
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are substantially below the current QCD sum rule results. 

Certain ratios of branching fractions can shed more light on the underlying hadronic dy- 
namics. The hadronic uncertainties on the two ratios 



vltlnZp+p] = °- 80 -o:K:i (ex P : 0.79 ± 0.17), (27) 



^11{bZ P / p \ = O-^KSS (exp: 0.046 ± 0.016) (28) 

are determined almost entirely by spectator scattering such that a larger ratio implies that 
this mechanism is more important. On the other hand, 



BrAv(5 -► p-u) 
BrAv(B p-p [ 



1.2 ^"""(O) 



4^(0) 



2 



(o.68±8:gg±8:SS) (ex p: o.^il) m 



provides insight on the ratio of the B — > oj to B — ► p form factor, as indicated by the above 
dependence on this ratio. More interesting information could be obtained from the pp final 
states, once the semi-leptonic B — > plv spectrum is measured more accurately near q 2 = 0. 

The predicted direct CP asymmetries are either very uncertain (colour-suppressed modes), 
often preferring only one or the other sign of the asymmetry, or rather small (colour-allowed 
decays). The small asymmetries for the colour-allowed decays follow from the dominance of 
the longitudinal polarisation amplitude combined with the smallness of the penguin amplitude. 
The available measurements are consistent with small or vanishing asymmetries, but do not 
allow to draw further conclusions at this moment. 



5.2. Longitudinal amplitudes and the determination of ol (7) from Sl 

For phenomenological studies it is often convenient to parameterise the decay amplitudes by 
hadronic amplitudes that can be directly determined from data. In the limit of isospin sym- 
metry and neglecting electroweak penguin contributions, the B — ► pp amplitude system is 
conventionally written in terms of complex graphical "tree", "colour-suppressed tree" and 
"penguin" amplitudes, 

^A-^p-p" = (T h + C h )e- l \ 

A B^p +P - =T h e-^ + P\ (30) 
- A B^pOpO = ~ pk - 

A similar set of equations applies to the B — > tttt system. This amounts to five real hadronic pa- 
rameters per helicity amplitude. Given 7 they can be extracted from the three helicity-specific 
branching fractions, the direct CP asymmetry in B° — ► (p + p~)h, and the time-dependent CP 
asymmetry S^ p , and compared to theoretical calculations. 

We calculate these quantities in QCD factorisation, where the main contributions to T, C and 
P come from the coefficients a±, «2 and a^+P^, respectively. In the following discussion we will 
only consider the longitudinal amplitudes, drop the helicity index and write C = \T\ x rce lS ° , 
P = \T\ x rpe lSp . The results are given in Table 4, which also compares the pp to the tttt 
system. The errors in this Table are from hadronic parameters, re, 5c, Sp do not depend 
on CKM parameters. The uncertainty from \V u b\ and V c b can be included noting that T is 
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B — > pp B — > tttt 



|T|/(KT 8 GcV" 1 ) 




o i c+0.58 
z - io -0.55 


r c 


o.3oi8:i* 


o-57±8:2S 


S C 


(-ioll)° 


(-411)° 


r P 


12+ 002 


0.421^ 


5 P 


(8l^ 4 )° 


(0lg)° 



Table 4: Amplitude parameters of the /jp and tttt systems. Errors do not include CKM param- 
eter uncertainties. 



proportional to \ V u b\, rp to |V^,|/| V u b\. Similar results have been presented in [48]. Numerical 
differences arise from a different choice of input parameters (for instance, here we use the same 
value QAe tlpA = 0.6e~ l40 ° for pions and p mesons) and the inclusion of spectator-scattering 
effects at next-to-next-to- leading order in [48]. The values reported here and in [48] provide a 
good description of all available tttt and pp observables within uncertainties of the calculation 
with the exception of the direct CP asymmetry in B° — > tt + tt~ , which is predicted to be smaller 
than what is observed (by the BELLE experiment). 

It is interesting to understand the difference between the pp and tttt system. The value of 
|T|, which controls the absolute magnitude of the colour- allowed branching fractions, is larger 
for p mesons, because the product f p A^^ p {Q) of decay constant and form factor is larger than 
for pions, see Section 4. The second important difference is the smaller penguin-to-tree ratio rp 
for p mesons, which follows from the absence of the power suppressed but "chirally-enhanced" 
scalar penguin amplitude. This also causes rc to differ, because T = a.\ + ol\ + . . ., and 
C = Q2 — ol\ + . . ., and explains why the branching fraction of B° — ► p + p~ is about four times 
larger than B° — ► tt + tt~, while those of B° — > p°p° and B° — > 7r°7r° are about equal. 

We now turn to the determination of 7 (or a) from time-dependent CP violation. The two 
CP asymmetries are defined through 

TUBy r py-)-TdBy r pyn = 

r L (BV(t)^p+p-) + r L (B°(t)^p+p-) ; ' y ' 

where Am > is the mass difference of the two neutral B meson mass eigenstates. We obtain 

Cf = 0.0271^™! (exp: -0.11 ±0.13), (32) 

in good agreement with experiment. As emphasized in [25, 31], the asymmetries Sf are par- 
ticularly suited to determine the CKM phase in the framework of QCD factorisation, because 
hadronic uncertainty enters only in the penguin-correction term to Sf, and the dependence on 
the strong phase 5p comes through cos 5p in very good approximation. Thus, like in no other 
observable, the hadronic uncertainty is much smaller than the dependence on the CKM phase 
7. This is especially true for the pp system, where rp is small, so that 

S P L P = sin 2a + 2 r P cos S P sin 7 cos 2a + 0{rp) (33) 

with a = 7r — (3 — 7. 

Figure 2 shows S^ p as a function of 7 in the range of interest. The band quantifies the 
theoretical uncertainty. From the intersection of this band with the measured value —0.06 ± 
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30 40 50 60 70 80 90 100 
y/degrees 



Figure 2: Theoretical result for the longitudinal time-dependent CP-asymmetry parameter S^ p 
as a function of the CKM angle 7. The horizontal band indicates the current exper- 
imental constraint [44, 49]. 

0.18 [44, 49], and given (3 = (21.2 ± 1.0)° [50, 51, 52], we obtain 

7= (73.2±™)° or a = (85.6j£f)°, (34) 

where the theoretical error alone is only ±3°. The value of 7 obtained in this way is remarkably 
consistent with the one from QCD factorisation calculations of the S'-parameters of the 7r + 7r~ 
and 7r ± p =F final states [25, 53], and presently provides the most accurate direct determination 
of 7. It is also consistent with [54], where instead of a theoretical calculation of P/T one 
uses SU(3) symmetry to relate the penguin amplitude to the longitudinal branching fraction 
of B~ —* p~K*°, and with other determinations of 7 (a) from B — > pp decays [44, 49]. 

5.3. Polarisation observables 

We now study the transverse-helicity contributions, which manifest themselves in polarisa- 
tion observables. As explained before, model-dependent effects such as non-factorisation of 
spectator-scattering and penguin annihilation can make a strong impact on transverse ampli- 
tudes, and therefore our results suffer from larger uncertainties. Here, this specifically concerns 
the five colour-suppressed modes. On the other hand, as polarisation observables like fi involve 
ratios, other uncertainties are often reduced. Specifically, CKM factors often drop out approx- 
imately, and form factors only enter in form of ratios (for example transverse/longitudinal), 
when only one form factor contribution is present or strongly dominant in an amplitude. 

Our results for the longitudinal polarisation fraction and the corresponding CP asymmetry 
are shown in Table 5. As expected the colour-allowed tree-dominated decay modes are pre- 
dicted to have fi near 1 with errors in the (5 — 10)% range. Their longitudinal CP asymmetries 
are predicted not to exceed 10%. Again the situation is very different for the colour-suppressed 
decays. With the exception of B° — > cap there is still a preference for significant longitudi- 
nal polarisation, but much smaller values can be obtained within theoretical errors. The large 
downward uncertainty is entirely due to the model-dependence in the spectator-scattering con- 
tribution to the negative helicity amplitude (parameterised by Xh in the QCD factorisation 
formalism). The longitudinal CP asymmetries are also very uncertain. 
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f L I percent 



Theory 



Aqp I percent 



Experiment Theory 
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Table 5: Longitudinal polarisation fraction and the corresponding CP asymmetry for tree- 
dominated decays. Experimental values are taken from [42, 43, 44, 45, 46, 47]. 




Figure 3: Correlation of branching fraction and longitudinal polarisation fraction in B — > p p , 
illustrated using 5000 randomly chosen points in our parameter space. The shaded 
area corresponds to the BABAR measurement [45], the highlighted point to our 
default values. 
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4>\\ I degrees A</>y / degrees (<j>^ — <j>±) / degrees 

B- - P- P ° -5±8±i^ ±2 
+l±|±i? +4±1±1 ±0 

B° — > p°/9° no prediction no prediction ±3 

b-^u P - — 9±|±i| +18±«±if ±3 
50 -29±t±H -45;^ ±1 

— > no prediction no prediction ±3 

B a ^K*+p- -l±i±|i +7™ ±0 

2? s — ► K*°p° no prediction no prediction ±2 
2? s — > K*°a; no prediction no prediction ±3 

Table 6: Predictions for other polarisation observables 

The theoretical predictions of Jl are compatible with the present experimental data where 
available. One notices, however, that the pattern of 1 — /l for the two colour-allowed pp final 
states seen by experiment appears to be opposite to the theoretical one though perhaps not 
significantly. We therefore calculate 

r + o = \Z f f L ^p\ = o-46±8:82±S:S («p= 2.75^). (35) 

The theoretical upper limit of r + o ~ 2 is attained when spectator-scattering is minimal 
(Xh = — 1), in which case /l{p p°) ~ 0.90 and /l(p + P~) ~ 0.95. The branching fraction 
and longitudinal polarisation fraction of the colour-suppressed decay B° — > p°p° are theo- 
retically allowed to lie within large ranges. To investigate the question whether there exist 
(theoretical) correlations between the two observables, we perform a random scan through the 
theory parameter space. The result is displayed in Figure 3. It shows that while there is a 
preference for smaller branching fractions than in our default prediction, there is no obvious 
correlation between the branching fraction and Jl- 

The remaining six polarisation observables can be taken to be /y — f±, Aq P — A^ p , and the 
phase differences between the transverse helicity amplitudes and the longitudinal amplitude. 
As explained in Section 3, f±, Aq P , <fij_, A<p± are expected to be approximately equal to /m, 

A^p, 0||, A0||. We find indeed that /y — f± is always below 2%. The last column of Table 6 
quantifies the expectation of equal 4>\\ and 4>±: the difference of the two does not exceed a few 
degrees. We should point out, however, that this calculation relies on the assumption that 
the positive-helicity is not substantially different from its magnitude in naive factorisation. In 
view of this the smallness of (f>^ — (fr± should be interpreted as the statement that no concrete 
dynamical mechanism is known that could produce a larger difference. The phase observables 
0|| and the corresponding CP asymmetry A^y are shown in the second and third column of 
Table 6. For most of the colour-suppressed decays the theoretical uncertainty is above 270°, in 
which case we conclude that no useful theoretical prediction is possible. For the colour-allowed 
modes, we obtain reasonably accurate results, which could be compared to experiment, once 
a complete angular analysis is performed. 
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Figure 4: Transverse-to- longitudinal amplitude ratios for B° — ► ET*(/>. The contours indicate 
the dependence of the negative-helicity theory prediction on annihilation (left) and 
spectator scattering (right) parameters, with all other input fixed at central val- 
ues. Additionally, the values implied by current measurements [26, 27, 29] for both 
helicities are shown. 



6. Penguin-dominated decays 

The 14 decay modes that we study in this section are characterised by the dominant role of the 
colour-allowed QCD penguin amplitude a 1 ^ = + /3f, which includes a penguin-annihilation 
term. Of these the 11 AS = 1 modes have branching fractions up to 10 5 , and some of them 
have already been studied extensively experimentally including polarisation. 

Due to their common dominant amplitude the theoretical errors in this class of decays are 
common to all representatives. As explained in Section 3.2, the negative-helicity penguin 
amplitude is particularly uncertain due to a potentially large penguin weak annihilation 
contribution [6]. In addition, non-factorisation of spectator scattering also affects the transverse 
amplitude of final states containing u or <j> mesons, mostly through the flavour-singlet penguin 
amplitude af - . An important issue of the subsequent analysis will be whether theoretical 
calculations are compatible with the observation of large transverse polarisation, and whether 
uncertainties can be controlled to the point that useful predictions can be made. 

6.1. The B — ► 4>K* system and the transverse penguin amplitude 

We begin with a discussion of the B — ► <pK* modes. A complete angular analysis is available for 
B — ► (j)K*° [26, 29], which allows us to extract the complex amplitude ratios A±/Ao from data. 
This is shown in Figure 4, which compares this result to the theoretical calculation of A-/Ao- 
(The experimental result for A+/Ao is in very good agreement with the expectation that the 
plus-helicity amplitude should be strongly suppressed.) The left plot in the figure shows the 
theoretical range from a variation of the uncertainties in weak annihilation alone (parameter 
Xa), the right plot displays the same information for spectator scattering (parameter Xh). 
Since all values for inside the contour are theoretically allowed for A—/Aq, it is evident that 
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Table 7: Comparison of theoretical results for observables from the full angular analysis of 
B~ -> cj)K*- and 5° -> (f)K*° with experimental results from [26, 27, 29]. 



theory does not require the amplitude ratio to be small. While it does not make accurate 
predictions, it is natural that penguin-dominated decays exhibit large transverse polarisation. 
This is confirmed by comparing the first column of numbers in Table 7 with the measurements 
in the fourth column. We find very good agreement of our results with data but with very large 
uncertainties. We also note that all observables related to the positive-helicity amplitude (the 
difference of || and ± observables) are predicted to be very small. So are the CP asymmetries, 

(s) 

since the doubly CKM-suppressed amplitude proportional to A„ ; does not exceed a few percent. 
Unless experiments find unexpectedly large values for any of these observables, the interesting 
ones are the branching fraction, fi and the phase (/>». 

We now explore a strategy where the variation of input parameters or the transverse penguin 
amplitude oP^ is constrained by data in order to improve the predictions for other observables 
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Figure 5: Predicted branching fraction and longitudinal polarisation fraction of the pure pen- 
guin decay B~ — > K*°p~ (left and center) and longitudinal polarisation fraction of 
B~ — > K*~p° (right) as a function of the parameter qa of our annihilation model, 
showing only the theory uncertainty from variation of the phase <fA- The horizontal 
band represents current experimental values [50]. The lines within the theory band 
indicate our default choices for qa and tpA- 

and decay modes. We assume that this amplitude is approximately the same for all decay 
modes, in accordance with factorisation calculations. Figure 4 allows explanations for the 
large negative- helicity amplitude based on spectator scattering or weak annihilation (or both). 
We favour the second option, since an enhancement of spectator scattering would lead to large 
transverse polarisation for the colour-suppressed tree decays which is not observed for the p°p° 
final state, see Section 5. Further support for this option comes from the pK* final states, 
which do not involve a|, and which are therefore much less sensitive to spectator scattering. 
Figure 5 displays three observables from the pK* system as a function of the strength of 
weak annihilation, qa- It can be seen that it must be non-zero, with a favoured range around 
qa ~ 0.6. This is consistent with a fit of qa e llfA to the (j)K*° data, which suggests 

qa = 0.5 ± 0.2 exp . <p A = (-43±19 exp .)°, (36) 

excluding other theory uncertainties. This coincidence motivates the default input parameter 
choice qa e lipA = 0.6e~' 40 ° adopted in Section 4. The second column in Table 7 shows the 
theoretical prediction when the variation of these parameters is reduced to qa = 0.6 ± 0.2 and 
ip A = (—40 ±10)° as suggested by (36) and Figure 5. The central values of these results remain 
the same by construction, but the hadronic uncertainties are considerably reduced. 

Rather than relying on our model-dependent parameterisation of the weak annihilation am- 
plitude to fit data, we prefer the point of view that is theoretically unreliable and should 
be taken from data. Thus, instead of qa^ a we fit . Neglecting CP violation in the 
(pK* system in accordance with the present data and theoretical expectations, the transverse 
amplitude can be expressed as 

A- = A^X^P^, (37) 

where, neglecting small coefficients, P_ « a1~(K*(j)) + a^(K*(j)). Calculating A- from the 
data with the overall phase adapted so that the phases of the longitudinal helicity amplitudes 
from theory and data match, we obtain 

P5*^ = (-0.081 ± 0.002(exp)±°;°°S(th)) + i (0.026 ± 0.011(exp)i^|(th)), (38) 

(s) 

where the theoretical error includes the uncertainties from A^^^Xc only. Since we already 
assumed that a large enhancement of is not a favoured option, we obtain 6l\ by subtracting 
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B-^K*-(f) lO.llg;!^ 2 10.4+°|i|; 2 9.7 ±1.5 

B0^K*°4> 9.3$i±l¥ Q - Q -oi-fl 9.5 ±0.8 

B-^K*-u 2.4^i 2: | 2.3i^t^ <3.4 

B°^K*°u 2.0±°}±fi 1.9±g-i±J-| <4.2 

B-^K*°p- 5.9lg-^-? 5.8±g"i±f-S 9.2 ±1.5 

B -^K*-p° 4.5±i;g±?;2 4.5±i;g±i;g < 6.1 

R0 . + c c+l.7+5.7 c 4 +1.7+2.6 ^10 

a ^ IX p °-°-l.5-2.9 °- 4 -1.5-1.5 < iZ 

B°^K*°p° 2.4i°; 2 l^ 2.3±g;?±S;l 5.6 ±1.6 

B.^K-K*+ 9.1±H±™ 8.0l 2 ;^|I n/a 

£> s -> 9.1+0.6+11.3 7.9+0-4+4.3 n/a 

2i-8±i:i±??:S w-B^g^ 1 i4.o;? : ° 

B-^K*°K*- 0.5±g-?±g-| 0.5i°- 2 t°- 2 <71 

B0^K*0K*° 0.6±°;}i°;f 0.6±8:i±g;l < 22 

g L=1 £^ 0.4±g;j±g;l O.algjjlg;? n/a 

Table 8: CP-averaged branching fractions of AS 1 = 1 and AD = 1 penguin-dominated 2? — > V 
decays. Experimental values are taken from [26, 27, 29, 42, 55, 56, 57]. 



the default value of . Rounding numbers, this results in 

a c ~ = (-0.08 ± 0.02) + i (0.03 ± 0.02). (39) 

This is the input to the theoretical predictions shown in the third column of Table 7 and 
the columns labeled u a^~ from data" in later tables in this section. This procedure provides 
another considerable reduction of hadronic uncertainties. 

6.2. Branching fractions, direct CP asymmetries and polarisation 

We now discuss the complete set of final states, where the a c A amplitude is dominant. The CP- 
averaged branching fractions, direct CP asymmetries and longitudinal polarisation fractions 
are given in Tables 8, 9 and 10, respectively. The two columns of numbers in these tables 
represent the result with our default inputs with uncertainties, and the result using (39) as 
input. In general, we regard the second result as our "best" prediction. However, one should 
be aware that it depends on the assumption of final-state independence of a%, and the present 
experimental data. 

The two uiK* modes are predicted to have 4-5 times smaller branching fractions than the 
(j)K* modes, which is consistent with experimental upper limits. Following the notation of the 
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Table 9: Direct CP asymmetries of AS* = 1 and AD = 1 penguin-dominated B 
Experimental values are taken from [26, 27, 29, 42, 55]. 



VV decays. 



appendix of [25] and suppressing helicity labels h, the decay amplitudes read 
V2A B -^uK*- = ^-K'uj I 2a 3 + 5 pu 02] + A uK , [a p 4 + 5 pu Qi] , 



V2A 



A K*co I 2 «3 + 6 PU «2] + A wR » 



(40) 



4' 



where we do not show numerically irrelevant amplitudes. The smaller branching fraction is 
simply a consequence of a relative factor of y/2 in the amplitude, and the smaller B — > w form 
factors multiplying the dominant coefficient (l^^-T^M^. | 2 ~ 2.5). As both factors cancel 
out in the polarisation fraction, Jl is predicted to be similar for tuK* and (f>K*. However, 
contrary to the <pK* system, the presence of tree amplitudes on 2 in (40) allows for sizeable 
CP asymmetries as can be seen from Table 9. 

The four K*p final states are the VV equivalents to the much discussed ttK final states. 
Their amplitudes are given by 
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(41) 
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Table 10: Longitudinal polarisation fraction of AS* = 1 and AD = 1 penguin-dominated B 
VV decays. Experimental values are taken from [26, 27, 29, 42, 55, 56]. 



They are particularly interesting, because the colour-allowed electroweak penguin amplitude 
"•Tew pl a y s an important role, in particular for polarisation. We discuss this point separately 
in Section 6.3. Here we note with respect to Table 8 that the branching fractions seem to 
be systematically below the measurements. Since the B~ — > K*°p~ mode is a pure penguin 
decay, proportional to the B — > p form factor, this is problematic. A larger B — > p form factor 
is not an option, since this would be in conflict with the observed B — > pp branching fractions. 
A larger value of d^' would also increase the <pK* branching fractions, unless the penguin 
amplitude is highly non-universal, or the B — > K* form factors are smaller than assumed, or 
one arranges a cancellation between a / 1 and a^ h in 4>K*. Rather than pursuing any of these 
options, we leave this issue as a potential problem for the QCD factorisation approach (or the 
input parameter set). The magnitude and sign of direct CP asymmetries is again related to 
the presence or absence of tree amplitudes that can interfere with the leading QCD penguin. 
Therefore, we predict negligible Acp for the pure-penguin mode K*°p~ , while asymmetries up 
to about 30% are possible for some of the other modes. 

The list of AS* = 1 penguin-dominated decays terminates with three B s modes. Here the 
flavour topology allows the penguin-annihilation amplitude which turns out to be the first 
subdominant contribution besides a.% for the longitudinal amplitude, but is negligible for the 
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transverse ones. Including these coefficients, the simplified amplitude expressions read 

(42) 

•ABs-^K*- K*+ = -^-K*K*04 + ^K*K* [<^4 + + $pu a l] j 



"4 + «3 - \ a lEW + ^4 



The relative factor of two in the 00 amplitude compared to the others leads to a particularly 
large branching fraction for this decay, even though the enhancement is somewhat reduced due 
to the destructive interference of the electroweak penguin and colour-suppressed QCD penguin 
amplitude in the longitudinal amplitude, see Table 2. The interference is constructive for the 
negative-helicity amplitude, leading to smaller fi for B s — > 00 compared to the other two B s 
decays. Sizeable CP asymmetries are only expected for the decay B s — > K*~K* + with a tree 
contribution. 

We briefly examine the three AD = 1 penguin channels, which have small branching frac- 
tions. The first subdominant contribution in these decays comes from the 0% up-penguin 
amplitude. In theoretical calculations the strong phases of a\ and are correlated to a cer- 
tain extent. This information is lost for the transverse amplitude when a%~ is taken from data 
and a%~ is calculated. This explains why the error on the CP asymmetries increases for these 
decays from the first to the second columns of numbers in Table 9. 

We conclude this subsection by providing the remaining polarisation observables in Table 11. 
As we see from Figure 4, the transverse phase observable 0y is very sensitive on penguin weak 
annihilation, and therefore uncertain, but it is expected to be nearly the same for all 14 decay 
modes. This is clearly also true when we fit ol\ to data, but in this case the errors are 
much smaller. For the polarisation CP asymmetries, the same qualitative statements as for 
the full direct CP asymmetries are valid, i. e. they can be significantly different from zero only 
in decays with a tree contribution or in the AD = 1 modes. Only rough estimates for these 
asymmetries can be given. We do not display the observables that vanish when the positive- 
helicity amplitude is zero. We find that |/m — f±\ < 4% for all modes, and that 0m — 0j_ 
and A0|| — A0^ never exceed ±2°. As discussed before, these statements should be taken 
with some caution, since the posit ive-helicity amplitude is estimated in the naive-factorisation 
approximation. 

6.3. The B — >■ pK* system and the electromagnetic penguin effect 

The pK* final states are particularly interesting for an investigation of electroweak penguin 
effects, since the suppression of the leading VV QCD penguin amplitude makes the electroweak 
penguin amplitude sizeable in comparison (about 50%). Moreover, as can be seen from (41), 
the electroweak penguin enters the amplitudes in three different combinations, 

A pK .a% h + k-A R . p la$* w , k = 1,0,-1, (43) 

allowing various kinds of interferences. Another interesting point is that the pattern of in- 
terference is opposite for the longitudinal and negative-helicity amplitude, since the sign of 
a 3 EW ^ s different from a^ w . This comes from an additional, power-enhanced contribution 
to «3'e W [20], such that 

a 3,EW|excl. a 3,EW ~ a 3,EW| excl . 3^ W 7 m 2 ' ^> 
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Table 11: Theoretical results for other polarisation observables in penguin-dominated B — > VV 
decays. 

which changes the real part from — O.OIOIq'qq^ to the value +0.015^q'qq3 given in Table 2. 
Since the term proportional to the Wilson coefficient of the electromagnetic dipole operator, 
Cj^, is the largest contribution to the negative- helicity electroweak penguin amplitude, the 
interference patterns (43) are sensitive to possible anomalous contributions to Cf^, including 
its phase. 

In Table 12 we compare selected observables for the two final states involving 03'^, when 
the extra term in (44) is excluded, to the default (included) and data. We note that already 
in the "excluded" results, the longitudinal polarisation fractions of the pK* final states are 
predicted to differ such that f L (K*-p°) > f L (K*°p-) > f L (K*°p°). This follows from the 
large longitudinal electroweak penguin contribution. The transverse electromagnetic dipole 
effect amplifies the hierarchy among the three fi predictions. The current experimental data 
confirm the first inequality, but the second is not seen. 

Similar to what has been done for the irK system in [25, 58, 59, 60, 61, 62], one can con- 
struct decay rate ratios that highlight the electroweak penguin contribution. However, in 
contrast to 7rK, one cannot expand in small amplitude ratios; the suppression of the VV 
QCD penguin amplitude makes these ratios too large. In [20] a few amplitude ratios related 
to the transverse polarisation decay rates have been discussed, and their dependence on the 
electromagnetic dipole operator has been emphasized. The default input of the present anal- 
ysis is similar to [20]. The amplitude ratio = P/f w /Ph used in [20] is approximately 
equal to 3a c 3 '^ w /(2a% h ). We now calculate Re(f£ w ) = — 0.25±?;1| [+0.15^;^] compared to 
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BrAv /10~ 6 


4.5 


5.4 


< 6.1 


2.4 


1.4 


5.6 ±1.6 


h/% 


84 


70 




22 


37 


57 ± 12 


^CP / % 


16 


14 


2011 


-15 


-24 


9 ±19 



Table 12: Predicted branching fraction, longitudinal polarisation and direct CP asymmetry of 
the two pK* final states sensitive to the electroweak penguin amplitude with the 
power-enhanced transverse contribution proportional to CV 7 included or excluded 
(theoretical errors in earlier tables). Experimental results for comparison (exp.). 



Default Ph from data 
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0.36 




2.34 




0.15 


5_ 


1.53 


S'_ 


0.51 


Si 


2.98 



Table 13: Longitudinal and negative-helicity CP-averaged decay rate ratios as defined in (45). 

Numbers in square brackets refer to the unrealistic scenario when the electromagnetic 
dipole effect is switched off. 

Re(p5 w ) = —0.23 ± 0.08 [+0.141q q 5 ] in [20], where the numbers in square brackets refer to 
the (unrealistic) scenario when the electromagnetic dipole effect is switched off. Only a small 
part of the large range of theoretically allowed values given in the first number is in fact com- 
patible with the observed branching fractions, so that the determination of the longitudinal 
and negative-helicity from theory alone is not the optimal approach. Therefore, in [20] 
the QCD penguin amplitudes Po,- have been obtained from the branching and longitudinal 
polarisation fraction of the pure penguin decay B~ — > K*°p~, and the assumption that the 
phase of does not exceed ±30° has been made. The experimental data that goes into this 
analysis has changed since; most importantly, the current smaller value of fL(K*°p~) implies a 
larger transverse QCD penguin amplitude. Repeating the fit in [20] with current data, we now 
obtain Re(p5 w ) = — 0.211qq| [+0.11±0.06], which reduces the impact of electroweak penguins 
by about 20% relative to the calculation, but dramatically improves on the theoretical error. 
We define the helicity-specific CP-averaged decay rate ratios [20] 

_ 2T h (p°K*°) , _ 2^0^*-) 

Sh T h {p-K«) ' Sh - T h (p-K*) ' (45) 

and S h = Sh/S' h . It should be emphasised that the CP-average of helicity-specific decay rates 
is not the same as the CP-average of polarisation fractions fh- When the standard variables are 
used, the relation involves CP asymmetries. The S*-observables defined above are better suited 
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to an investigation of helicity-specific effects. Experimentally they can be determined from the 
same data as the standard observables, thus avoiding unfolding complicated correlations in the 
errors of CP asymmetries, branching and polarisation fractions. In Table 13 we summarise 
our theoretical predictions for these ratios, providing results for both cases (default input and 
Ph from data), and with the electromagnetic dipole effect switched off for comparison. Note 
that all ratios would equal 1, if the QCD penguin amplitude was really dominant. The largely 
different numbers illustrate the impact of electroweak penguins in these decays, both for the 
longitudinal and transverse amplitudes. One also observes that the electromagnetic dipole 
operator contribution is essential in the transverse case. Using data to fit the leading QCD 
penguin amplitude from the pure penguin mode is obviously crucial to discriminate the effect. 
It is most pronounced in the S_ ratio, where the electroweak penguin amplitude enters the 
relevant amplitudes with opposite sign, see (41). 

7. Other B — ► VV decays 

In this section we briefly discuss the remaining 11 decay modes consisting of electroweak or 
QCD flavour-singlet penguin-dominated decays and decays that can only proceed via weak 
annihilation. None of these modes depends on the QCD penguin amplitude d^, so we give 
the numerical results according to the general procedure defined in Section 4 and used for the 
analysis of the tree-dominated modes in Section 5. 

7.1. Other penguin-dominated decays 

These five decays modes are characterised by an interplay of the flavour-singlet QCD penguin 
amplitude a§, the colour-allowed electroweak penguin amplitude 03 ew> an< ^ the c °l° ur " an d 
CKM-suppressed tree amplitude 02- All three amplitudes are significantly smaller than a\, 
hence these decays have small branching fractions compared to the AS = 1 penguin-dominated 
decays. 

We first consider the three AD = 1 B~, B° decay modes in this class, whose amplitudes, 
following the notation of the appendix of [25] and neglecting the helicity labels, are given by 





— A p ^ 


p 

«3 — 


1 P 1 

2 a 3,EW 




= A p( p 


P 

a 3 ~ 


1 P - 

2 a 3,EW 




= 


P 

«3 " 


1 P 1 

2 a 3,EW 



The three decay modes have an identical amplitude structure. It can be seen from Table 2 
that the QCD penguin and electroweak penguin contributions are of similar magnitude for 
the dominant longitudinal amplitude, and interfer constructively. Nevertheless, the branching 
fraction of the first decay is expected to be only (1 — 3) • 10~ 8 , and the second and third are 
a factor of two and three smaller, respectively. The negative helicity amplitude is dominated 
by the electromagnetic penguin contribution due to the power-enhanced effect discussed in 
[20] and Section 3.2.3, but it is not larger than the longitudinal amplitude. The longitudinal 
polarisation fraction fi is predicted in the range 0.7...1. All CP-violating observables are 
small, since QCD factorisation calculations do not produce a significant strong phase between 
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fr, / percent 
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B s - #° 
-B s — ► 


n 4n +0. 12+0.25 
u - 4U -0.10-0.04 
n 1f| +0.05+0.48 
— u.uo — v.iz 


01+3+9 
8i -4-12 

no prediction 


177+6+9 

L < '-7-15 

iqi+1+51 
101 i inn 

— 1 — 1UU 




Aqp j percent 


^cp / P ercen t 


A0M / degrees 


5 S - #° 
-Bs — > cj0 


iq+5+56 

O+3+102 
°-3-56 


i 1+4+10 
11 -3 -8 

no prediction 


i c+4+19 
i0 -4-14 

1 O+5+80 

J-J-4-72 



Table 14: Predictions for S s — > w0 and -B s — > </>/?. 



the two terms with different weak phases, a% — EW and a% — EW . Our results for the 
branching fractions are compatible with those of [63]. 

The two B s modes in this class have larger branching fractions, since they are governed by 
AS* = 1 transitions. The amplitudes are given by (helicity labels omitted) 



$pu «2 + ^ a 3,EW 

5 pu a 2 + 2o| + ^«3 EW 



(47) 



Due to a partial cancellation between the QCD and electroweak penguin contributions, the 
CKM-suppressed tree amplitude a 2 is the largest partial amplitude in the second decay. The 
roles of tree and penguin amplitudes are reversed in the first decay. Our results for the various 
observables are summarised in Table 14. They are rather uncertain for B s — ► uxf), where 
the non-factorisation of transverse spectator-scattering is important, such that often no useful 
prediction can be obtained. It is also worth noting here that we have assumed ideal mixing 
throughout this paper, such that the u meson has no ss component. Since the amplitude for 
B s — > <j)(j) is an order of magnitude larger than for B s — > uj(f> even a small mixing angle of about 
5° could make a significant difference in the results for B s — > tocf). The polarisation observables 
of B s — > p o are determined by the power-enhanced contributions from the electromagnetic 
dipole operator, which dominates the transverse electroweak penguin amplitude. Similar to 
B — ► /xfT* decays [20], this contribution changes the sign of af EW relative to naive factorisation, 
and hence changes 0y by almost 180°. 

7.2. Pure weak annihilation decays 

Branching fraction estimates for the six decay modes that can proceed only via weak annihila- 
tion are given in Table 15. Since QCD factorisation does not provide a solid prediction of the 
annihilation amplitudes, these numbers should be regarded as estimates within the adopted 
annihilation model. Measurements of these decay modes would result in useful checks of this 
model. We do not present other observables for the annihilation modes, because the calcula- 
tions are too crude to provide quantitative results. The following qualitative conclusions can, 
however, be drawn: the decay amplitudes of these modes depend only on the tree annihilation 
amplitude (3\, and the penguin annihilation amplitude 0% (and a corresponding electroweak 
penguin annihilation amplitude), but not on 0^. The calculation of the transverse annihilation 
amplitudes shows an enhancement only for as discussed in Section 3.2.4, while all others 
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rO k*~K*+ n nq+ - 05 + - 12 
21 u - uy -0.03-0.10 

B° ^<M> < 0.03 

R __ n Q A +0-03 +0.60 

B s ^ p p 0.34_ 03 _ 38 

R y rp rp 1 7 +0.01 +0.30 

D s ^pp U.l '_o.01-0.19 

£? s -► _>p° < 0.01 

d __, , ,, , n i 1+0.01+0.20 

_f s -> U - iJ --0.01-0.12 

Table 15: Branching fraction estimates for pure weak-annihilation modes. 



respect the hierarchy (10) of helicity amplitudes. Thus we expect /_ to be in the range 0.8 ... 1 
and f± w /|| for all modes listed in Table 15. Experimental tests of these expectations would 
be very interesting for the understanding of annihilation dynamics, but require rather large 
data samples. 



8. Conclusion 

In this paper we performed a comprehensive analysis of the 34 B decays to two vector mesons 
(and their CP conjugates). Together with [25], where the PP and PV final states were 
discussed, this completes the phenomenology of two-body decays in QCD factorisation in 
next-to-leading order. In comparison with the PP and PV final states, the VV ones are much 
more uncertain. This is due to a potentially large negative-helicity penguin weak-annihilation 
amplitude pointed out in [6], but also due to the non-factorisation of spectator-scattering for 
the transverse amplitudes, which has a particularly large effect on colour-suppressed partial 
amplitudes. Our main results are summarised as follows. 

Results related to tree- dominated decays 

We obtain a very good description of the pp system including the p°p° final state. As a general 
rule, all colour-suppressed tree decays are poorly predicted due to the non-factorisation of 
transverse spectator-scattering. The observed large longitudinal polarisation /_(p p°) suggests, 
however, that this effect is not as large as it could be. In contrast, the colour-allowed tree- 
dominated decays are the theoretically best predicted B — > VV modes, and should all show 
/_ near 1. We also find a small longitudinal QCD penguin amplitude, which makes the time- 
dependent CP asymmetry S p £ an ideal observable to determine the CKM angle 7. We find 

7 =(73.2i™)° or a = (85.6^)°, (48) 

where the theoretical error alone is only ±3°. 

Results related to penguin- dominated decays 

The penguin-dominated decays are plagued by the weak-annihilation uncertainty. While it is 
natural to obtain an equal amount of transverse and longitudinal polarisation, many observ- 
ables can only be predicted if at least some information is taken from data. In our analysis, we 
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explored the possibility of replacing the calculated negative-helicity penguin amplitude by a fit 
from (j)K* data. While weak annihilation remains the most plausible dynamical explanation 
for significant transverse polarisation, other options exist: spectator-scattering may enhance 
the transverse flavour-singlet QCD penguin amplitude, and the electroweak penguin ampli- 
tude receives a large contribution from the electromagnetic dipole operator [20]. The latter 
effect is expected to be most prominent in the pK* system, and we propose to measure certain 
helicity-specific decay rate ratios to isolate it. 

Finally, we comment on the possibility to uncover the helicity structure of the weak interac- 
tions through polarisation studies in B — > VV decays. For instance, the presence of tensor 
operators changes the helicity-amplitude hierarchy (10). According to the theoretical picture 
that emerges from our study, this appears to be very difficult as far as the hierarchy between 
the negative-helicity and longitudinal amplitude is concerned, since it is already weak or vi- 
olated in the QCD penguin amplitudes due to Standard Model QCD dynamics. The tree 
amplitudes cannot receive large new contributions, since similar effects should then be seen 
in semi-leptonic decays. This leaves the colour- allowed electroweak penguin amplitude, which 
is theoretically well-controlled, providing further motivation for the investigation of the pK* 
system. Regarding the hierarchy between the negative- and positive-helcity amplitudes, there 
is currently no indication of its violation, neither from the (j)K* data nor from theory. How- 
ever, one should be aware that the factorisation properties of the positive-helicity amplitude 
are virtually unknown. 
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Note added 

After publication of the preprint version of this paper, we have been made aware of talks by 
A. Kagan, where a fit to the B — > <fiK* system similar to the one discussed in Section 6.1 was 
reported, arriving at similar conlcusions on the magnitude of transverse weak annihilation. 

A. Appendix 

In this appendix we collect the formulae for the hard-scattering functions relevant to the 
negative-helicity amplitudes. 

A.l. Light-cone projection 

Up to twist-3, six two-particle light-cone distribution amplitudes are relevant for vector mesons. 
We neglect three-particle qqg amplitudes. Consistency then requires that one adopts the 
Wandzura-Wilczek relations among the two-particle amplitudes, reducing the input to two 
amplitudes <E>y = <j>\\ and 4>± [64]. 
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The two-particle light-cone projection operator on the longitudinal polarisation state of the 
vector meson is given as given in [25], Section 2.3. To obtain the projector on the transverse 

polarisation states in the helicity basis, we use the result for M\_ from [33], and insert e± = e T 
to obtain 
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if v m v 
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where p is the momentum of the vector meson, and we assigned the parton momenta 
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(49) 
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to the quark and antiquark constituents, respectively. The collinear approximation l± = 
must be applied only after the projection. Using the Wandzura-Wilczek relations [64] and 
defining [6] 



J u 



dv 
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o v 



(51) 



the transverse helicity projectors simplify to 



4 

ifvmv J 



V 



+ <Pb{u) 



7"(1 =F7 5 ) + ujf(l T75) 



(52) 



7"(1 ±75) - ^^(1 ±75) 



l ± =0 



Many previous calculations of transverse polarisation amplitudes [4, 5, 12, 16] used an incorrect 
expression for the projector, where the transverse-derivative terms are dropped. The authors 
of [13] used (49) without employing the Wandzura-Wilczek relations. The projector is not 
given explicitly in [6]. 

A.2. The af h coefficients 

The leading order results for the coefficients af h , needed to calculate amplitudes according 
to (15), (16), reproduce the results from naive factorisation. At next-to-leading order they 
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receive contributions from one-loop vertex corrections, penguin and dipole operator insertion 
topologies and hard spectator interaction terms. We assemble these in the form 



a, 



?;±i 
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Cj±\ Cpa s 



N?(V 2 ) 

4tt 2 

v?(y 2 ) + — HfrvM) 



+ P l h ' P (V 2 ), 



(53) 



where the upper (lower) signs apply when i is odd (even), and it is understood that the 
superscript 'p' is to be ommited for i G {1,2}. 

The results for h = correspond to those given in [25] for PV final states with obvious 
replacements of P by V, so we will only give explicit results for the negative-helicity amplitude. 
The leading-order coefficient, corresponding to naive factorisation, is simply 



N~(V 2 



0, ie{6,8} 

1, else. 



(54) 



As is well-known, there is no leading-order contribution from (S — P) <S> (S + P) operators for 
vector mesons. 

The negative-helicity vertex corrections read 



V~(V 2 



f 

Jo 



dy 4>b2(y) 



dy 4>a2(y) 



121n^-18+<7 T (y) 

-121n-^+ 6-g T (y) 
A 1 



0. 



i £ {1,2,3,4,9,10} 

ie{5,7} 

% e {6,8} 



with (y = 1 — y; similarly for other convolution variables below) 
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The function gxiy) differs from the corresponding function g(y) in the longitudinal amplitude 
only by a single term. 

Corrections from penguin contractions and dipole-operator insertions are present for i G 
{4,7-10} at order a s (a em for the electroweak penguin amplitudes i = 7-10). We find 
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P~' P (V 2 ) = P7' P (V 2 ) = 0, and 
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(60) 



where G(s,y) is the usual penguin function (see, for instance, [31]). Notice that unlike the 
longitudinal case, there is no contribution from the dipole operators Qs g and Qjy to P^' P (V 2 ), 
P^ ,p (y 2 ), for which we confirm the result given in [6]. The two terms P^~' P (V 2 ), P 9 ~' P (V 2 ) exhibit 
an enhancement by (m b /AQcr>) 2 proportional to Cfjjf due to the factor mBm b /m 2 , which alters 
the naive power-counting (10) to (21). See also [20]. We neglect the small contributions from 
the electroweak penguin operators Q7-10 to the penguin coefficients. 
The spectator-scattering contribution is given by 
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for i e {1,2,3,4,9,10}, 
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for i £ {6, 8}. Our expressions for the spectator-scattering kernels are simpler than previously 
published results [4, 5, 12, 13, 16], and differ even when the Wandzura-Wilzcek relations are 
used to simplify those results. The B meson light-cone distribution amplitude enters the 
kernels via the parameter Xb defined in [23]. The factors in front of the integrals are of 
order 1 in the Aqcd/uh, counting except for i = 6,8, where there is an extra ms/ Aqcb- 
enhancement. Thus, the hard-spectator scattering contribution to a p G g is formally leading over 
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the naive-factorisation amplitude, although not numerically, because Hq g are multiplied by 
small Wilson coefficients. Since a? g contribute to a decay amplitude in the product r 



since rY 2 ~ 



v 2 p— 

X fl 6,8' 



and 



^QCD/^b according to (17), we conclude that the suppression of the (pseudo) scalar 
penguin amplitudes relative to the V— A ones, a^ - , is absent for the negative helicity amplitude. 

Another qualitative difference to the longitudinal amplitude is that the x-integrals in (61), 
(62) are divergent, since the integrand is too singular near x = 1. Thus factorisation for the 
negative-helicity amplitude breaks down even at leading order in the heavy-quark expansion 
due to the non-factorisation of spectator-scattering. To estimate this contribution, we extract 
the logarithmic divergence by applying a "plus-prescription" to the integrand, and replace 
the large logarithm by a phenomenological parameter [25]. To estimate the endpoint 
behaviour, we note that the asymptotic distribution amplitudes are (j)^~{x) — > 6xx, <j) a {x) — ► 
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finite 



As specified in (25), we use a simple model where we treat Xjj as an unknown complex param- 
eter universal to all Hi(ViV 2 ) with magnitude around ln(m.B/AQCD) because the logarithmic 
infrared divergence has its origin in a soft gluon interaction with the spectator quark and can 
therefore be expected to be regulated at a physical scale of order Aqcd- 

A. 3. Weak annihilation contributions (b? ,h coefficients) 

Weak- annihilation is parameterised by a set of amplitudes b^(V\V 2 ). The leading contributions 
can be assembled from a few basic building blocks as shown in [25, 31]. The corresponding 
formulae also hold for the helicity amplitudes in B — > VV decays, hence we only summarise 
these building blocks here. 

For the longitudinal case h = 0, only a few signs change in comparison with the known 
results for B — > PP or PV . We find 
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f f 

and A L f = = 0. The non-vanishing transverse building blocks are 
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A^ = ira s 



dxdy 



2m 



A 



f- 



4>ai(x)<p2(y) _ ^2 r vi <t>j{x)4>b2{y) 

m 2 x xy(l — xy) mi x xy(l — xy) 
2mi V2 <t)ai{x)4>^{y) 2m 2 Vl {x)<j> h2 (y) 



f 1 {' 

na s J dxdy I /\-'— + /■; 



m 2 



xy^ 



x 2 y 



(66) 



i+ 



2m±m 2 



m\ 



dxdy { 4> a i(x)<j) a2 (y) 



xy 3 (1 — xy) 2 y 2 {l — xy) 



A' + = 7ra s 2mi ™ 2 J fady ^ 0/i[ _(,r)o h2 (ii) 



m. 



x 



x 



x 3 y (1 — xy) 2 x 2 (l — xy) 



Again, the convolution integrals exhibit logarithmic and even linear infrared divergences, which 
we extract into unknown complex quantities using the prescriptions 



1 du 



u 



Xa, 



f 

Jo 



du 



Inu 



u 



1 du 
u 2 



X L . 



(67) 



As for Xh, we assume both Xa and Xl to be universal to all VV final states with magnitudes 
around ln(m;>/AQCD) and toj,/Aqcd, respectively. 

Since all non- vanishing building blocks contain such divergences, making the treatment of an- 
nihilation rather model dependent, we further simplify our results by evaluating the convolution 
integrals with asymptotic distribution amplitudes <£y(u) = (f>\_{u) = 6uu, 4> a {u) = 4>b(u) = 3m 2 , 
& v (u) = 3(u — u) and obtain the expressions 



4-° 



4° « 18vra 



X A -4 + yJ + r^r^{X A -2) 

4° « 18ttq s (r^ + r^) {-X 2 A + 2X A - 4 + , 
4° « 18vra s {r^ - r^) {2X A - 1) (2 - A A ) 



(68) 
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for the nonvanishing longitudinal contributions, and 

4+ « 4+ « 18™ s ^ (2X1 - 3X A + 6 - |vr 2 ) , 

.\ mimo / 1 „ 5 7r 2 \ 

B V 7 (69) 

4" « 18™. f I^r? - ^rA (X\ - 2X A + 2) , 

4- « 18^a s f ^rj + (2X1 ~ $X A + 3) 

for the transverse ones. As expected, all annihilation contributions are suppressed by at least 
one power of Aqcd/^-b compared to the form factor contributions, with 4'°) 4° an( ^ 4 + 2 
carrying an additional explicit suppression factor. Where a comparison is possible, these results 
agree with [6]. 
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